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Abstract 

We present the Lagrangian whose corresponding action is the trace K action for General 
Relativity. Although this Lagrangian is second order in the derivatives, it has no second 
order time derivatives and its behavior at space infinity in the asymptotically flat case is 
identical to other alternative Lagrangians for General Relativity, like the gamma-gamma 
Lagrangian used by Einstein. We develop some elements of the variational principle for field 
theories with boundaries, and apply them to second order Lagrangians, where we stablish 
the conditions — proposition 1 — for the conservation of the Noether charges. From this 
general approach a pre-symplectic form is naturally obtained that features two terms, one 
from the bulk and another from the boundary. When applied to the trace K Lagrangian, 
we recover a pre-symplectic form first introduced using a different approach. We prove that 
all diffeomorphisms satisfying certain restrictions at the boundary — that keep room for a 
realization of the Poincare group — will yield Noether conserved charges. In particular, the 
computation of the total energy gives, in the asymptotically flat case, the ADM result. 



1 Introduction 

Either way, boundary conditions from boundary terms or boundary terms for boundary con- 
ditions; boundary terms for the action and boundary conditions for the field configurations go 
together. Whereas from a purely logical point of view the boundary conditions emerge as a con- 
sequence of the boundary terms, from the technical side the situation is often the opposite: one 
must look for the boundary terms that are necessary to implement the desired — or acceptable — 
boundary conditions. 

This discussion applies directly to General Relativity (GR). It is well known that the bound- 
ary conditions required for the correct application of the variational principle in GR depend 
upon the boundary terms exhibited by the action. Since divergence terms do not alter Einstein 
equations of motion, there is some freedom to write down Lagrangians for GR that differ in 
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some divergence terms from the original Einstein-Hilbert proposal. These divergence terms lead 
to different boundary terms for the action. To be specific, consider the original Einstein-Hilbert 
action (integration is in a 4-volume of spacetime) 



Seh = / d^x J\g\R . (1) 

In this case, the boundary conditions for the application of the variational principle involve the 
fixation of some derivatives of the metric at the boundary. Instead, subtraction of a divergence 
from (HI) allows to write the "gamma-gamma" action, first used by Einstein, 

Sr = A y^5^-(r^,r^, - r^^r^j , (2) 

which is first order in the spacetime derivatives and that has the possible advantage that the 
variational principle only requires to fix the metric at the boundary. 

In field theory, the association -by means of the Noether theorem- of symmetries at the level 
of the variational principle with conserved currents makes the role of boundary terms (leading 
to boundary conditions for the fields), either in the Lagrangian or the Hamiltonian formulation, 
essential for the conservation of charges, as shown in the pioneering work of |jl| for GR (see also 
recent work in |Q, ^, Q). The role of boundary terms has also been stressed in recent years by 
the introduction of the concept of quasilocal charges in GR (see Q for general references, see 
also I, g, 1). 

Many actions, all differing from Einstein-Hilbert's in boundary terms, have been used in the 
literature. With no aim of completeness, let us mention contributions in 10, 11, 12, 14, 



15, 16, 17, 19, 20| and many references therein. Among them, and using the terminology 
of g, the "trace K" action |ll|, g offers some singular features, because of its geometric 
character, that deserve full interest. This action is roughly written as 



Sk = / A J\9\R-'^ / A J\-f\K , (3) 

V Jg^^ V 

where 8%^ is the boundary of the 4-volume 7 is the determinant of the 3-metric induced 
on the boundary and K is its extrinsic curvature. When the boundary is not smooth, and 
discontinuities arise in the vector field orthonormal to the boundary, the second term in (^) 
must be understood ||l^ as including delta-like contributions from the "joints" between the 
smooth elements of the boundary. The detailed form of 5"^: is given in section 3. 

The purpose of the present paper is to analyse the Lagrangian leading to the action (^) and 
to provide with the theoretical framework for actions of this type. In section 2 we introduce 
some notation and useful formulas. The Lagrangian for the trace K action is obtained in section 
3. In section 4 we develop the formal theory for the variational principle and the Noether 
symmetries for second order Lagrangians in field theories with boundaries. This theory is applied 
in section 5 to the trace K Lagrangian and we compute the total energy in the particular case of 
asymptotically flat spaces. Conclusions are presented in section 6 and the appendix is devoted 
to prove some expressions used in the text. 
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2 Some notation and formulas 

Here we introduce some notation and formulas to be used in the next sections. We consider that 
our spacetime coordinates correspond to a standard 3 + 1 decomposition, including in particular 
that 5°° < and g'' > 0, that is, no surface = constant is tangent to the light cone. Working 
with this type of coordinates means no restriction, at least infinitesimally, on the gauge freedom, 
because any infinitesimal diffcomorphism preserves these conditions for the metric components. 
The determinant of the 4-mctric g^i, will be denoted by g. 

Given the 3-surface = constant, its orthonormal vector n(//) is defined by the relation 

and by requiring that it is pointing towards increasing values of the coordinate x^. The coefficient 
^{/j.) is just a sign, ^(fj,) = 77^^^, where 7)^1, is the Minkowski metric with positive signature 
(-,+,+,+). Then, 



9^' 



^A9^\ 

Let us point out that n(/Li) is not a true vector field for it fails to transform as such under 
the diffeomorphisms that do not preserve the foliation defined by the 3-surfaces x'^ = constant. 
We compute in the appendix the deviation from the vector behavior of the transformation of 
n(/x) under diffeomorphisms. 

The 3-metric induced at the 3-surface = constant is gab, with a,b = 0,1,2,3 except /i. 
Its inverse matrix is given by 7"^(At); the components of this matrix are the non identically 
vanishing components of 

j^if,) := - = gf>^ - e(/x)n^(/x)n-(M). (4) 

The determinant det gab (for a,b = 0,1,2,3 except fj,) will be denoted 7(/x). There is the 
relationship 

55'^'^=7(/^)- 

Consider, for u ^ fx, the 2-surface x^ = constant, x^ = constant. The induced metric on it 
is gAB, with A,B = 1, 2, 3 except /x, v. Its inverse metric will be written as 7'^^(/xz/), and the 
determinant det gAB ='■ One can show: 

gig^^g'''' - {g^^) = gg^^Y'^ii^) = 7(m^^)- 

Note also that 

Y'if^) = g'-^'ii - = g'-'ii - im{^)q\H\ (5) 

where we have defined the scalar products 
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The extrinsic curvature for the surface x'^ = constant is given by 
which is equivalent to 



where >Cn(^) is the Lie derivative under n(/i). 

The trace of the extrinsic curvature K{fj,) := j"'^{ij,)Kah{fi) may be written as 



Kin) = i{^i)C^(^^){ln'^\g\) = (,ip)n'f^{n) , 

where n"jj{^j) is the covariant derivative of n(/x), V^n'^ifJ'), with the Riemannian connexion, as 
if n(;u) were a true vector. 

3 The Lagrangian for the trace K action 



In this section we use techniques introduced in ||21[ (see alsop2|, |23||). Let us first take a look 
on the boundary conditions imposed by the variational principle for the Einstein-Hilbert La- 
grangian. A general variation for Ceh gives 



6CEH = -G''''Sg^, + ^\g\g''''6R^,, (6) 
where R^^, stands for the Ricci tensor and 

G^" := vS(i?^^ - ^g'^'R), 

with R the scalar curvature. 

As it is well known, the last term in (0) is a divergence. It can be written as 



gig'-'SR,,, = d,,{gr5V%) , (7) 

where 5^^^ := y/\g\{gP'' - g^^'^'^d^^). The next step is to write Q g^^^ST'^u as (we continue 
with the convention of indices a,b = 0, 1, 2, 3 except fj,) 

gr^r;, = -gab6P''\ii) + dA^\g^^5{^)) , (8) 



where Pabil^) (for a given /i, notice that we raise and lower indices with gab,! (fJ-)) is a gener- 
alisation to any /x of the ADM |24| momenta for /i = 0: 

PabifJ-) = \l\-f{lj)\{gabK{ll) - KabifJ-)). 

Therefore, the divergence term in (^ is 

d.igr^K'^) = -d,{gabif^)5P'^\f^)) + d,d,(^^\g'^'^5i^)^ 

= -d,{gabif^)6P^''{f^)) + ld,d,{^\ [g^^6{^) + <7^^5(^)) ). (9) 
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Finally, the term acted upon by the derivatives df^di, has the equivalent expression: 

-\f\9\ (g^'^H—) + g'^'Si—)) = (\f\9\\/Cifj-)C{i^)gf'f'g''A 6 \ , , , = 



/|7(HI'5a(H , (10) 

where the angles a{iiv) are defined as 

Q(Oi) = arcsinh(g(Oi)) , Ci{ij) = aicsm.{q{ij)) , 

and we have used 



V(i-c(/^)ew(^^)) 

Therefore 

5Seh = d^x 6£eh = d^x [-G^'^Sg^, - d^{gab{fi)6P^Hl^)) + df,d,{^\^ifiu)\6aifii^)) 

(11) 

Consider the volume %^ being a 4-cube, having as boundary elements 3-surfaces defined by 
the constancy of one of the coordinates x^. Then the second term in the integrand will give 
a contribution from the faces and the third a contribution from the joints, these joints being 
defined by the constancy of two coordinates. In order for Seh to be a differentiable functional 
1 25 1, and its variation to consequently lead to Einstein equations, we need to control the fields 
and the variations at the boundary in such a way that 

SSeh = ^ = 0. 

As we have already said, this control at the boundary involves some derivatives of the 
metric components. Now we have the specifics: the variation of P°'^{fj,) must vanish at the 
faces x^ = constant, and the variation of a{fiiy) must vanish at the joints between the faces 
x^ = constant and x'^ = constant. The physical meaning of this restrictions on the variations 
is unclear. 

Instead, a simple addition of a divergence term to Ceh will give a more reasonable control of 
the variations at the boundary. Taking into account that (7af,(/i)P"''(/i) = 2^/\J{fl)\K{fl), define 
the new Lagrangian (K is for the extrinsic curvature) 



Ck := Ceh + 25^(y^^K(/i)) - d^d^i^/hii^Hf^u)) . (12) 
The action (^ is, by definition, Sk '■= Jiy d'h: Ck- Then, 

5Sk = j^^ A 5Ck = A (^-Gf^'dgf,, + df,{P'''{^i)6gM) - ^^^,{a{^iu)6^Jh{^^\)^ 



(13) 

So for Sk to be a differentiable functional we must require, besides the customary vanishing 
of the variations of the 3- metric induced on the initial and final equal-time 3-surfaces, the 
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supplementary condition of the vanishing of the variations of the 3-metric induced on any spatial 
3-face of the boundary. These requirements already guarantee the fixation of the variations of 
the determinant ^{^^v) at the joints. In fact, to be more precise, the vanishing of the variations of 
the 3-metric is only required for finite boundaries; in the case of asymptotically fiat spaces, where 
we let the spatial elements of the boundary go to the space infinity (r — > oo), we must control 
the r — > oo behavior of the fields and its allowed variations so that there is no contribution to 
5Sk- More on this later. 

Notice the difference with the boundary conditions that one finds in the case of the gamma- 
gamma Lagrangian, where the vanishing of the variations of the 4-metric is required. 

Expression (12) is the Lagrangian for the action (P). We notice that the contributions from 



the joints are all included. The role of these contributions has been explained in [15|. Early 



computations can be found in [22 



Two features of (|T^) are worth being mentioned immediately. First, as it happens generally 
with the proposals to modify Ceh through divergence terms, Ck is not a truly scalar density. 
This fact does not opposes to its physical applicability, as it is argued in a parallel context in (2^ 
for the gamma-gamma Lagrangian. Ck gives indeed the trace K action, with all its geometric 
meaning, but for coordinates adapted to the boundary -or viceversa: for boundaries adapted to 
the coordinates-, in such a way that the elements of the boundary correspond to the constancy 
of the value of some coordinate. A typical boundary may be the 4-cube, already used, but a 
cylinder whose top and bottom faces are equal-time 3-surfaces, and whose lateral face is defined 
by the constancy of a single radial coordinate, is a well adapted boundary as well. In the cylinder 
case the {ij) is for space indices) contributions from the last term in ( p^ ) disappear in Sk 
because "the boundary of a boundary is zero" . 

The second observation is that, unlike the "gamma-gamma" Lagrangian (Q), Ck is not a first 
order Lagrangian. This is a point not sufficiently recog nised in the literature 0, ||, |3, H]. It 
is proved in the appendix that the second order contributions to Ck are as follows: 



= (quadratic terms in the first derivatives of the metric) — a{iii')d^di,\J\^{nv)\, (14) 
but note a key difference with the Einstein-Hilbert Lagrangian: Ck has no second order time 



derivatives (the sum fi,u in (14) only contributes for ^jl ^ v). In this sense Ck has an intermediate 
place between Ceh (with second order time derivatives) and Cy (with only first order spacetime 
derivatives). In the next section we will introduce some notation and results for theories with 
second order Lagrangians to focus later on Lagrangians of the type of Ck- The advantage 
of not having second order time derivatives will become clear when we analyse the boundary 
conditions. Also, we will see in the asymptotically flat case that the long distance behavior of 
Ck improves crucially that of Ceh- 

4 Second order Lagrangians: variational principle for field the- 
ories with boundaries 



Here we will consider a generic second order Lagrangian density function C with dependences: 

C{(j),4)f„(i)f,y) 
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(p denotes the whole set of fields (a new index could be introduced but it will be unnecessary). 
(pfj, stands for d^cp, and := d^dycj). If, as it happens with Ck, only terms with ^ ^ v appear 
in the second derivatives, some of the expressions simplify somewhat. 
The Lagrangian functional is 



for some spatial 3-volume, and the action functional is 

S[4>] = [ dt L[(l),4>,4>] = [ d^x C{(t),(l)i„4>f,u), 
with '^V = [to,ti] X The functional differentiation of a general S is given p^] by 



5S= / (15) 



(m is a condensed notation for any number of partial derivatives with respect to any spacetime 
coordinate), that must be compared with 

6S = J A 6C = J A ^S^(-^ (16) 

in order to define the functional derivatives j^§^- We will ignore, in the way of computing 
variations with respect to the second order derivatives, the symmetry (p^u = (p^^. This means 
that we compute independently, for instance, and 5^^- It is useful then to define the 

"complete" derivative — that takes into account the truly independent variables — , 

d^S 5^S dS 5S 



uip^u •J'i'u^i uyj^^ ui^y^ 

for /i 7^ I/, and 



S^S 5S 



'-'tip. "V^/i/i 



for the rest. 

Comparing ([l5|) and (|l^) gives, 

5S dC dC dC 

d(f> d(p dcpfj. d<p^iy 

^ = (19) 

The first line displays the Euler-Lagrange derivatives, yielding the equations of motion when set 
to zero. 

Analogous variations can be computed for L[(j), (p, (p], taking into account that (p, (p and (p are 
independent arguments for L. It is convenient to define 

dC dC dC 5S 5^S 5L 5^L 

+ 



d(pi ' dcpQi d(pio' 
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For the sake of completeness we give the relationship between functional derivatives for S and 
L. They are 

6S 6L ^ ,SL. - - ,SL, 5S 5L ,^L. 5S 5L _ ,5L. 
TT = 7T-^o(--)+5o5o(-y), -y = - 29o(-y), TTT = TTT " ^0 ) > 



and 



5S 6L 



6S 



6L 



Now we will explore the requirement of differentiability coming from the variational principle. 
Using (HI) 



SS 



V 



V 



5(t) 



+ 



dui 



6S 



93 V 



3y 
t=ti 

t=to 



t=tl 



+ 



-I t=to 



tl 



dt I dui 



to 



+ / A didj{ 

Uy 



6S 



+ 



3y 0000 



t=ti 
t=to 



(20) 



where d^V is the -spatial- boundary of and dai are the two-forms induced at the boundary 
by the application of Stokes theorem. 

The second term, the fourth term, and a piece within the sixth term, in the last equality, van- 
ish because 6(j)\tQ = 6(t>\t-i = as conditions imposed on the variations allowed by the variational 
principle. So in order for the variation of S to determine the Euler-Lagrange equations, 

5S 



6S = Q 



, 



we need 



ii r 
dt / dai 

to JdW 



5S 



6S 



+ 



d^x 



6S 



3V 



6<p) 



-i t=ti 



0. 



(21) 



t=to 



Note in the last term the presence of 6(p. This term complicates the setting of the boundary 



conditions for the fields. Formula (21) — or better, formula (po|), before the cancellations orig- 
inated from the variational principle — is responsible for the expression ( pT| ) obtained for the 
Einstein-Hilbert action. Things become a lot easier if we assume our Lagrangian not to contain 
second order time derivatives, which is the case for Ck- We will continue with this assumption. 



4.1 Second order Lagrangians with no second order time derivatives 

This implies that now L is L[(j),(j)] only. In this case, (p]) simplifies to 



dt 



d(Ti -r—ocp + dil-r— 



0, 



Ho Jd^V \<J(Pi 0(pij j 

and since the — finite — range for the time integration is arbitrary, we end up with 



0. 



(22) 



(23) 
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Condition ( |23D must be fulfilled, at any time, by the fields and its variations in order to 
comply with the variational principle. Two restrictions originate from (|2J 



A restriction on the space of field configurations — spatial boundary conditions for the 
fields. 

Consequently, a restriction on the allowed variations so that the space of field configurations 
is preserved under such variations. 



Both restrictions altogether must imply (|23|). Obviously the simplest restrictions one can think 
of are to fix the values of the fields at the spatial boundary so that its variations at the boundary 
vanish (maybe not all the fields need to be fixed at the spatial boundary, because some of the 
functional derivatives of S in ( p3| ) may vanish identically). This trivially complies with (p^). 



The restrictions set on the space of field configurations are subject to a consistency test: 
they must be compatible with the equations of motion. 

This compatibility is generally nontrivial, and may cause the rejection of some boundary con- 
ditions when they do not agree with the dynamics ( |2^ ) To implement the compatibility of 
(p3| ) with (^^, that must hold at any time, we can run, for a fixed time that can be taken as 
the initial time — when the initial conditions are set — , a Dirac-like algorithm of stabilisation 
of constraints. Then we can end up, in principle, with secondary boundary conditions, tertiary, 
etc., that play the role of gauge fixing constraints. We refer to pO[ | for a form of the standard 
Dirac algorithm for obtaining of the secondary, tertiary, etc., constraints, that makes no mention 
to tangency conditions. The idea that boundary conditions must be treated as Dirac constraints 



is clear from the developments above and has been put forward in |3f]. 

We have all things ready for the study of the implementation of the Noether symmetries for 
Ck- Before finishing this section, let us mention that the first line of equation (^) — or the 



^Sk rjah. 



more general version (|15| ) — , together with the variation of Ck in (P~3D, dictates the results: 

P'^\^,) (25) 



for a, 6 = 0, 1, 2, 3 except /j., and 

6Sk SSk 



-a(Hy'|7(H|7^^(H (26) 



for A, B = 0, 1, 2, 3 except fi, v. Note that (p^) is in agreement with (|T4|). These results will be 
used in section 5. 



^We leave aside the case when the boundary is an interface between two different physical regimes. 
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4.2 Noether currents and charges 

We first continue developing the general theory with the same assumptions of the preceding 
section: L is (p] and C is C{(l), (j), (j)i,<j)i, 4>ij) 

Now consider the Noether case. Suppose that our Lagrangian satisfies, under a certain 
variation 5d (not necessarily complying with (]2^)), 



5dC = a^F^ , (27) 
for certain functions F^. From the first line of (|20|), 



we get 

^^5d(^ = (f^^ - ^5d<I^ - d.{^6D<t>)\ , (28) 
which identifies the current 

J"" ■■=F^ -^^D4>-d,(^^bB4>)^d,A^\ (29) 

as an on shell conserved object, 

d,J^ = 0. (30) 

(on shell) 

Note that we have included in the unavoidable ambiguity of the addition of the divergence 
of an arbitrary antisymmetric object A^^ . When searching for a conserved charge, it will prove 
useful to take into account this ambiguity. Let us integrate ( |30D along the spatial 3-volume, 

9o / d\ f d(Ti J* = 0, (31) 

73 y Jd3v (on shell) 



SO the charge 



will be conserved on shell if 



It is advantageous to take, in J'^, 



Q := / d\ J° (32) 

J3y 



dOi J' = 0. 

Q3V (on shell) 



with B^^ another arbitrary antisymmetric object. Then the components of the current are 

j' = F'-^-^Jn<t^-dADct>-B^^) 

dtp 0(pi 

and 

d(pi d(pij 
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Thus, 

Jd^v Jd^v \ o<pi d(pij J 

and now (|2^ ) comes into play: the addition of the second and third terms in the right side 
vanish if the Eeld configurations and 6d comply with (23). Here we realize the relevance of 
the boundary conditions when considering the conservation of charge. Form now on we assume 
that the variational principle has restricted our space of field configurations and the allowed 
variations 6d defined on it in such way that ( p3|) is satisfied. Then the charge conservation (|32|) 
has the condition 



d'iV (on shell) 



whereas the charge itself takes the form 



Q:= f d\ J^= I d\ (f'- ^-^5d(^ - d^{^5D<t> - B^') 
J^v J'-'v \ 6(p S4>i J 



The presence of the object B^^ may be endowed to the ambiguity of in (pTD, under additions 
of the divergence of an arbitrary antisymmetric object. Absorbing dyB'^^ within and diB^'^ 
within F^ , we arrive at the following 

Proposition 1 Given a second order Lagrangian with no second order time derivatives. If 

1. the space of field configurations — including spatial boundary conditions for the fields — , 
together with the allowed variations defined on it, complies with (2S), 

2. the boundary conditions imposed on the field configurations are compatible with the equa- 
tions of motion ( 24 ), 

3. there is an allowed variation 6d such that 5dC = dfj_F^^ for some F^ , 
then, the charge 

Q = I d\ (FO - ^-SdcI^) - I da, {^Sd<P) (33) 
5(j) J a w 5(j)i 

is conserved if and only if the condition 



dai F' = , (34) 

93 y (on shell) 

holds. 

For infinite boundaries, as in the case of the limit r ^ 00 for asymptotically flat spaces, ( p^ ) 
is understood as 



hm / dai F' = . (35) 

'^^'^Jd^V (on shell) 

Note the generic form of the charge — conserved if the conditions of proposition 1 hold — in 
|): a bulk term and a boundary term, the existence of this last term having its origin in the 
dependence of the Lagrangian on the second order derivatives. 
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It is worth to remark that the Noether current conservation (^) stemming from (28) is a 
local property only of the equations of motion, and it is totally impervious to the modifications 
of the action by boundary terms. The Noether current is local, and for its definition ( p9[ ) it is 
completely irrelevant whether the variation 5d producing (p7|) is an allowed variation — in the 



sense of (|23D — or not. Let the action be Seh, Sr or Sk — each leading to different boundary 
conditions — , the current conservation equation (|30| ) will always be the same (but remember: 
there is an ambiguity in J^, to wit, the addition of a divergence of an arbitrary antisymmetric 
object). What then is the relevance of using one action or another? the answer is that the 
boundary terms in the action are indeed relevant to set the conditions ([2^). Selecting the action 
entails a selection of the space of field configurations and the allowed variations on it — though 
maybe not in a unique way. It is obvious then, looking at the proposition above, that the charge 
conservation crucially depends on the selected action. 

If the current conservation (^) holds but the conditions of Proposition 1 are not satisfied, 
then the equation ( |3l| ) will express the typical conservation equation that balances the rate of 
change in time of the total charge with the fiux of current through the boundary. 

4.3 Energy and (pre)symplectic structure 

Here we continue with a second order Lagrangian £. with no second order time derivatives. Since 
there is no explicit dependence of C on the coordinates, 6d(P = <pSt, with 6t an infinitesimal 
constant, is a Noether symmetry with F'^ = 6QC6t. Since F"^ = 0, the conservation of ( [3^ ) will 
always hold as long as conditions 1) and 2) in Proposition 1 are satisfied. With this proviso, the 



conserved quantity (33) is (up to a sign) the energy: 



E = f d\ (^<^ -C)+ j da, (^0) , (36) 
Jw 66 JdW 6 (hi 



that can also be expressed in terms of the Lagrangian functional, 

E= [ d\&^)+ [ da.U^^)-L[ 



Expression (3(:) is the formula for the energy for second order Lagrangians with no second 
order time derivatives. It generalises the common expression E = p^q^ — L{q, q), corresponding 
to the Legendre transformation in mechanics, where 

dL 
dq^ 



is the pullback to tangent space of the momentum in cotangent space (phase space) through the 
Legendre map. Some interesting consequences may be drawn from this generalisation. Summa- 
tion for an index k in mechanics becomes in field theory integration for the space coordinates 
plus summation for all the fields. Using this mechanical analogy, the role of the pullback Pkdq^ 
of the Liouville one-form is now played by 
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which can indeed be cahed the pullback of the Liouville form for this type of field theory. The 
pullback dpk A dq'' to tangent space of the symplectic two-form in phase space now becomes 

n= d^x [5^ ^6(P)+ / dai {5^ A 6^) . (37) 
J^v 54> Ja^v 6(t)i 

This structure can be symplectic (closed and maximal rank) or presymplectic (non maximal 
rank) depending on the presence of gauge symmetries in the theory. 

Expression (^^ also means that each canonical momentum has now two components: the 
bulk component and the boundary component. Their pullbacks to tangent space are, respec- 
tively, 

6S , , 

P(buik) = TT (38) 

and 

j bS 

P(boundary) = ■ (39) 

5 Application to C}q\ asymptotically flat spaces 



Here we follow Faddeev approach |2^. Asymptotically flat spacetimes correspond to physical 
situations where the gravitating masses and matter fields at finite times are effectively concen- 
trated in a finite region of space. Our spacetime will be a topologically simple manifold whose 
points can be parametrised by a system of four coordinates x^, — oo < < cxd, such that, in 
the limit r — > oo (r^ = {x}^'^ + (x^)^ + (x^)^) for finite time t = x^, the metric components 
satisfy the asymptotic conditions 

gtiu = Vtii^ + 0{^), dag^u = 0{^), 

= 0(1), ■■■,d(-^g,^ = 0{-^). (40) 



Conditions (40) amount to a partial gauge fixing, for the only acceptable changes of coor- 
dinates will be from now on those that preserve (|40|). To this end, for an infinitesimal change 
— > x''^ = x^ — e'^(x), we take to behave as 

e'^ = Lof'^x" + a'' + 0(-), 

r 

d,e^ = < + 0(^) 

9^^^'' = o(^),«>o 

d^"^^^' = O(^),«>0. (41) 

oj^^j^ is an infinitesimal Lorentz parameter, r]p^u:% + r]y^uj^p = 0, and a'^ is an infinitesimal 
translation of the coordinates. 
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Note that, under an infinitesimal diffeomorphism (^T]), 



= 0{r) X + ((r?,, + + 0{^^)) 

= 0{\), (42) 

in agreement with the asymptotic behavior (^). Also Sod^^^ g^iu = ^( ^.(J+i) )• 

Let us check that (pO| ) and (^) imply the boundary conditions (^). The integration in (23) 
is, for Ck, 

B.C. := I daJ^-^5^ + dAcP)] = f daJ'-^6^ + ld/-ll^ScP)] 
Jd^v \ 0(f)i d(pij J Jd'iv \ O0j 2 6(j)ij ) 

f 

where we have used (^) and (^6|). 

Since, when r — > oo, the area of the boundary d^V grows as O(r^), the asymptotic behavior 
of dH) is 

lim B.C. = limf da, {p''\i)6ga, - dMim^Hi3)\ " a{ij)d,{6^ \^m) 



da, {P''\i)5gab - \dM{ij)\IVm\l^''m9AB)^ 
dai [p^\i)5gab-dMij)^^\im)) > (43) 



;im 0(r2) (0(1) X 0(1) - 0(1) x 0{\) - 0{^) x 0(1; 
lim O(-) = , (44) 

r— +00 J- 



in agreement with ([2^). It remains now to check that there is consistency of the conditions (40) 



with the equations of motion ( p4| ) — for finite times. This can be verified using the Lagrangian 
equations of motion for Ck — or Ceh- The conditions of Proposition 1 are therefore satisfied. 

Let us apply to Ck the results of the preceding section. The bulk momentum is obtained 
from (p5[), 

^ = p^i(O) =: P^^ , (45) 

for i,i = 1, 2, 3; and the boundary momentum form (p6|), 

5Sk 



a(0i)^|7(0i)|7^''(0i) , (46) 



for A, i? = 1,2,3 except i. Let us find the Liouville and the presymplectic forms. The Liouville 
form is 



Pkdq^ = f A(-f<50)+ / da,{ 

J3v Sd) Jaw 6(pi 
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/3y Jd^v 
and the presymplectic form (p7[) is, accordingly. 



d\{P'^5gij) - / da, (a(007|7(0i)|7^''(0i)<55AB) 

- 2 / d^r^ (a(0i)5y^|7(0i)|) , (47) 



17 = / d^x{5 — ^ A (50) + / dai {6^^ A 
3y 5(h Jd^v 5(hi 



d^x (6P'^ A 6gii) -2 I da, (<5a(0i) A ^A/hfOi)!) . (48) 



This result was obtained in using the theory of symplectic relations 33, here it 

has been derived particular case of (p7|). 

According to the remarks produced at the end of subsection 4.2, since Ck differs from 
the scalar density Lagrangian Ceh by boundary terms, it is already guaranteed that Noether 
currents exist for all diffeomorphism symmetries. The point is that any Lagrangian £any differing 
from Ceh, a scalar density, by a divergence term, 

-Cany = CeH + d^,D>' , 

transforms, under an infinitesimal diffeomorphism 5d generated by the vector field e^O^ (not 
necessarily satisfying (pl|)), as 



5DC^ny = ^^,{e^'CEH + 5DD^'). 



5.1 The energy 



Let us consider now the energy for a theory described by Ck plus matter terms in an asymp- 
totically flat spacetime. We will assume that these matter terms are described by a first order 
scalar density matter Lagrangian that together with Ck define the total Lagrangian, and that 
the couplings with the metric are nonderivative. We will also assume that the matter terms 
satisfy the appropriate conditions at the boundary in order to comply with the variational prin- 
ciple. We will first work with pure gravity and at the end consider the changes introduced by 
the presence of matter. So we continue with the Lagrangian Ck- 
Expression ( |36| ) gives, for the energy density £ {E = Jsy d^x £) , 

£ = ^-l^-C + d,{^4>) 
Sep S(j)i 

= P'^gij -Ck- 2d,{a{Oi)do^Jh{Oi)\). (49) 



Introducing ( [l^ ) and recalling P^^gij = 2y/\^{0)\K{0), 

£ = P'^gij-^\R-2^,{^Hi)\K{i))-^o{P'^g^i, 



+ d,dj{^\^{ij)\a{i3)) + 2dido{^\^{Qi)\a{Qi)) - 2d,{a{Qi)d^^ \^{Qi)\) 



+ 5.5.(\/l7(u)l«(u')) + 25.(v/|7(Oi)|aoa(Oi)). (50) 
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Let us now rewrite some terms in (pO[). We show in the appendix, using the methods of |2 
that 

- P'^g,, + 2d,{^h{Oi)\doam = 2^\R^ - 25,(V^7°'^(0 (51) 

where Rq is a component of the Ricci tensor. On the other hand, using (p!0|), another piece in 
can be given a more convenient expression 

^lvm\dM^j) = \^\Ud,{^) + a"d,{^)\ 



\g\{T]^r'{i)+v]^r\j)) , (52) 

where f^g^" = has been used in the last equahty. Also, trivially, 



'm\K{i) = -^\g\^>^%i)V^^^. (53) 
All together. 



£ = ^\g\{2R', -R) + 2d,{^\g\^>^^{i)V%) - i^^\g\{r)^l^^ (i) + T]^^^\j)) 

+ a,(a(ii)a,-v^^)-2a,(v^7°^(i)n;.) 

= 2G^ + dM^j)dj^/\^\) + d,^^\{T)^r'i^-rl^^^ . (54) 



The first term in (|5j), a component of the Einstein tensor, vanishes on shell — it is a La- 
grangian constraint, part of the equations of motion. The total energy on shell is 

E = I d\ £\^^ gj^gii 



(on shell) J3y 



}}^l^^da, (a(u)a,v^^ + v^(rj^7^^(i)-r^:^7'^^(j))) . (55) 



Considering the asymptotic behavior (40), the contribution of the first term in (p^) vanishes 
because a{ij) = 0{^). The second term contributes 

E = lim / dai (T), - Vj-) = hm / dai {djgij - digjj) , (56) 

(on shell) r~^°°Jd3v ■' i~~^°°Jd3v 

which is the ADM |4| energy. 

The inclusion of matter, with the restrictions set at the beginning of this section, will change 
the term Gg to Gq — SttTq, where Tq is a component of the energy momentum tensor for the 
matter fields. Now Gq" - 8ttT(^ is a constraint that vanishes in virtue of Einstein equations. The 



asymptotic contribution in (56) remains the same. 



5.2 Other Noether charges for 

Space translations in the j direction (parameter in (p])) define for Ck the quantity F^j-^ in 
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that should satisfy p^ ) in Proposition 1. The r —>■ oo behavior of Ck can be easily deduced 
from the expression (p!^. The terms quadratic in the derivatives of the metric behave as 0{^) x 
0{^) = 0(4). The term a(/xi/)9^a,(/R^) has 



and 



Therefore 



a(^iu) ~ 0(-) 

r—*cx) f 



r^oo 



J — >oo 



lim / dai Ck = . (57) 
■ 'aw 



and so, 

This proves that the translational Noether symmetry leads to conserved quantities, the momenta. 

This result can be generalised: any diffeomorphism satisfying (41) has a conserved charge. 
To prove it we need to compute in (p7[) for an infinitesimal diffeomorphism 5d generated 



by e^dfj_. This computation is involved because Ck is not a scalar density. The definition (12) 
obviously gives 



5dCk = d^ie^'CEH) + 2d^ [5D{^\l{lj)\K{p))j - d^d, [5D{^\i{iJiiy)\a{ixiy))j . (58) 
but the behavior in the limit r ^ oo of the first term 

e^'CEH ^ 0(4) ' 

(with ~ Oir)) produces terms different from zero in (^). In fact this behavior is corrected 
by other terms in ( ^8|) but to single out these contributions is cumbersome. It is much more 
convenient to express Ck as the sum of the "gamma- gamma" Lagrangian Ct plus divergences, 
because the behavior of Cy in the limit r — > oo is much better than that of Ceh-, 

r— >oo 7"J 

So let us proceed to relate Ck with C^- Using (^3|), a trivial extension of (|^) to indices /uz/, 

5;.(v/|7(HI5.«(h) = -5^(v^(ri;,7'^^(A.) + r^.7^'^(z.))) 

= Id, (v^^K(^)) + (^(ri:,7^'^(M) - Kai^^^y))) m 

allows to write Ck in the equivalent form 

Ck = - d, (^a(Ha,v'|7(Hl) - (^^(r(;,7'^'^(/^) - Kai^^i^))) ■ m 

On the other hand, the "gamma-gamma" Lagrangian (^) can be written as 

Cr = MR-dJM{nag'''' -Kaf"")] , (61) 
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and therefore, using (^), 



-K 



(62) 



We continue with bo being an infinitesimal diffeomorphism generated by e^dn- Computing 



Sd^^t is a simple thing [26| because the deviation of Cy from the scalar density behavior is due 
exclusively to the two connexion coefficients present in the divergence term in (|6l|). Since the 
deviation of the connexion r(^^ from a tensorial behavior is an additive term e'^^, we can write 



5dCt = d. e^Cr 



and thus, for Ck, 



+ 6d 



(63) 



Now, from (p3), we deduce the object that needs to pass the test (p5| 



r = e'Cr-^J\g\g'"'e\^ + ^\g\g'''e^^,,-6D[a{^,.)^,^h{^,.)\ 

+ 5D(v^(ry(z)n'^(i)n'^(i)n.-r?(z.)n-(i.)nXi^)r^j) ■ (64) 

Now it is easy to check that each additive term in (|6^, under the conditions ( ^0|) and (|4l|), 
guarantees (35) and therefore 

lim / da^ F' = , 

thus proving that all diffeomorphisms satisfying ( |4l] ) give Noether conserved charges for asymp- 
totically flat spaces defined by the conditions (HO) . Note in particular the existence of conserved 



charges associated with the Poincare transformations included in (41). In this sense, (56) is 
just an example of the computation — the total energy — of one of the ten conserved charges 
corresponding to the Poincare group. 



6 Conclusions 

In this paper we have developed the formal theory for the variational principle and the Noether 
symmetries for second order Lagrangians in field theories with boundaries, with special emphasis 
in Lagrangians with no second order time derivatives. These developments lead to general 
expressions for some geometric objects, like the pullback to tangent space of the symplectic 
form in phase space. It is worth to remark that bulk terms and boundary terms contribute to 
these objects. Also the canonical momenta exhibit both a bulk piece and a boundary piece. 
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The Noether theorem is discussed for these theories. A clear distinction is made between 
conservation of currents, that is a local assertion only related to the equations of motion — but 
not to the possible boundary terms in the action — , and the conservation of charges, where the 
relation with the boundary terms appearing in the action is made transparent. The reason being 
that the boundary terms in the action are linked to the boundary conditions to be satisfied by 
the fields. These results are summarised in a proposition introduced in subsection 4.2. Let us 
mention also that the Noether charge exhibits a bulk piece and a boundary piece, the boundary 
piece owing its existence to the second order dependences in the Lagrangian. 

This framework is applied to the trace K Lagrangian for General Relativity. The trace 
K Lagrangian is obtained in section 3 and its applications to asymptotically flat spaces are 
studied in section 5. We observe that this Lagrangian is not first order in the derivatives but 
its asymptotic behavior is different from that of the Einstein-Hilbert Lagrangian and resembles 
that of first order Lagrangians like the gamma-gamma Lagrangian. We prove the conservation 
of charges for diffeomorphisms that preserve the boundary conditions for the metric tensor and, 
as a particular case, we obtain the ADM formula for the conserved energy. 

It is worth noting that we have based our approach in keeping at any moment the require- 
ments derived from the strict application of the variational principle, including in particular the 
differentiability of the action functional. These requirements lead to the imposition of boundary 
conditions on the field configurations. When we examine the conditions for the conservation of 
possible Noether charges the relevance of these boundary conditions, and hence of the variational 
principle, becomes transparent. 



Let us finally give a short review of the main results obtained in this paper. 
1. General Theory 

• In Proposition 1, section 4.2, we have given the necessary and sufficient conditions 
for the conservation of Noether charges for field theories with boundaries; these field 
theories being derived from action principles with second order Lagrangians that are 
free from second order time derivative terms. 



• We have shown that the Noether conserved charges, (|33|), exhibit a boundary piece 
that stems from the presence of the one-time-one-space derivative terms in the Lagran- 
gian. 

• We have derived a general formula, (^), for the (pre)symplectic form associated 
with this type of Lagrangians. This form exhibits a boundary piece that has the 
same origin as for the Noether charges. 

• Also, due to the presence of the one-time-one-space derivative terms in the La- 



grangian, we show that one is naturally led to define a bulk component, (38), as 
well as a boundary component, (^9[), for each of the momenta. 



2. Trace K Lagrangian 
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• We have neatly displayed, ([T^), the dependence of the Trace K Lagrangian on sec- 
ond order spacetime derivatives. We also show that the asymptotic behavior of this 
Lagrangian is similar to the one of the first order "gamma-gamma" Lagrangian, sub- 
stantially different from that of the Einstein-Hilbert Lagrangian. 

• We have applied our general formulas to compute the (pre)symplectic form associated 
with the Trace K Lagrangian, (^8|). Our results coincide with those obtained within 
the theory of symplectic relations. 

• We have given an expression for the energy density, (|5j), and we have computed the 
total energy for the asymptotically flat case, obtaining the ADM formula. 

• We have shown that, in the asymptotically flat case, all diffeomorphisms satisfying 
the boundary conditions -this includes the transformations that are asymptotically 
Poincare- give Noether conserved charges. 

7 Appendix 

7.1 Deviation of the transformation of n(/i) from the vector behavior 

Let 6d be the infinitesimal diffeomorphism generated by e^d^. Our task is to compute don^ilj) 
for the "vector" n(/i)whose components, n'^in), are defined through 



and we know 



Considering 



we obtain 



where we have introduced (5naive to symbolise a "naive" variation associated with vector behavior, 
that is, 

<5naiven"(/x) := c'' d^rf - n''{fi)d^e\ 
Next, to find S^n'^ilJ'), use 



and so. 



This expression can be arranged to give 



J.M. Pons, Boundary conditions and the trace K Lagrangian 



21 



which, using (^), is 

= 5naiven'^(/u)-n^(//)<5^7"'(/i)56e'^ , (65) 

for a,b = 0,1,2,3 except //. Equation (^) expresses the deviation of the transformation of 
n(^) from the "naive" vector behavior. Notice that this deviation differs from zero only for 
infinitesimal diffeomorphisms such that db€^ ^ 0; these are the diffeomorphisms that do not 
preserve the foliation of %^ in surfaces x'^ = constant. 



7.2 Proof of (|T4D 

Let us single out, for instance, the d^d^ terms in 



Ck = ^/\9\R + 2a^(/b(^K(/i)) - d^d,i^\j{fiu)\a{fiu)) . 
To this end, it is convenient to use the ADM decomposition 

^\R = Cadm - 2d,, (^^\{n\0)n>^^{0) - n^(0)nf;,(0))) 

where 

Cadm := -^hmnoiO)CR + KabiO)K^\0) - K''{0)) 

is free from dod^ terms, ^i? is the scalar curvature for the surface = constant. 

We will use the notation [something]^^ to isolate the additive terms in [something] that 
contain a derivative, and similarly for [something]^^^ . So 



[^i^]|o3 = -29o(v^[n^(0)nO,(0)-nO(0)nf,(0)]|3) 
- 2ds (V^[n^(0)nf,(0) - n\0)n^,iO)\^ 

+ 2^0 ([V^^i^(0)]|3) + 2ds ([./hmKi3%) - 2aoa3(v/|7(03)|a(03)) . (66) 



The first and third term in the right side cancel because n\{0) = and n\{0) = —K{0). 
On the other hand, the piece in the second term n'^(0)ni\(0) can be expressed as n^(0)n^^(0) = 
no{0)^^^ djn^ (0) , that has no do derivative. Therefore 



[Ck]io, = 2d3 (^[n3(0)nf,(0) + V^^i^(3) - ao(v/|7(03)ja(03))]|„ j , (67) 
Now, a little algebra gives 

[n'(0)nf,(0) + V^K(3)]|„ = ^^^nO(3)(/n^),o 



1 



vi + r(03) 



^|7(03)|aoa(03) , (68) 
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and, finally, plugging this result in 

[Ck]\os = -2a(03)ao93(Vl7(03)|) . (69) 

The ADM decomposition gives also an immediate proof that dodo terms are not present 
in £k- An extension of the ADM decomposition — which is associated with the n = 
coordinate — to any other coordinate helps to prove, along the same lines, that 

[^k]i,,=0, (70) 

and that 

[-Cii-] (second order terms) = -0:{|^l')^^,^„{^J \-f{nu)\) . (71) 



7.3 Proof of (|5TD 

Recalling (^), for 5 being the time derivative. 



= -9^,P'' + d.{^J\g\g''d,{^)) , (72) 



and recalling also ( |To[ ) (with 5 being the time derivative) 



Oi Oi 



we can write 



\9\ ( r5o(^) +5^^5o(|^) j = 2Vl7(00|5oa(Oi) , (73) 



Considering |21] the following identity with = (1,0,0,0), that originates from the defi- 
nition of the Riemann tensor, 

f^, = V^V.X'^-i?,^/X^ , (75) 

then 



= ^JJ\^\{\/^'XP -V^X^^)) +2J\^\R{X\ (76) 



and, for X^' = (1, 0, 0, 0) and p = 0, 



= di{^J\g\{S7'x'-V^X'))+2^\g\Ro'. (77) 

Then, 



Oi 



g,,P^-2d,i^h{Oi)\doam) = -^/^-of-^ - 9,(5-9o(|-)) 



-2V^< - a, Jl^m'X^ - V^X') + g''do{^)) 



2y^< + 25. (vS7°^(0n^) ■ (78) 
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